Quantum Nondemolition Squeezing of a Nanomechanical Resonator 
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We show that the nanoresonator position can be squeezed 
significantly below the ground state level by measuring the 
nanoresonator with a quantum point contact or a single- 
electron transistor and applying a periodic voltage across the 
detector. The mechanism of squeezing is basically a general- 
ization of quantum nondemolition measurement of an oscilla- 
tor to the case of continuous measurement by a weakly cou- 
pled detector. The quantum feedback is necessary to prevent 
the "heating" due to measurement back-action. We also dis- 
cuss a procedure of experimental verification of the squeezed 
state. 



I. INTRODUCTION 

Quantum nondemolition (QND) measurements [1~4] 
were proposed as a way to overcome the so-called stan- 
dard quantum limit for the measurement sensitivity, 
which arises due to quantum back-action of a detector. 
The general idea of a QND measurement is to avoid mea- 
suring (or obtaining any information on) the magnitude 
conjugated to the magnitude of interest, and therefore to 
avoid the corresponding back-action. An important im- 
plementation of this idea is the "stroboscopic" measure- 
ment of an oscillator position [2,3]. Suppose the position 
xi is measured (instantaneously) with a finite precision 
Aa;, which necessarily disturbs the momentum according 
to the Heisenberg uncertainty principle Ap > h/2Ax. 
Normally this momentum change would affect the result 
of the next position measurement X2 and would limit the 
accuracy for the position difference a;2 — xi , leading to the 
standard quantum limit for this magnitude. However, if 
the second measurement is performed exactly one oscil- 
lation period after the first one, the oscillator returns 
to its initial state, and therefore the momentum change 
does not affect the accuracy of X2 — xi measurement. 
Such stroboscopic measurement gives no information re- 
lated to the momentum, and this is exactly the reason 
why the effect of quantum back-action is avoided [1-4]. 

In terms of applications, the QND measurements have 
been mainly discussed in relation to measurement of very 
weak classical forces, in particular gravitational waves 
(see, e.g., [5,6]). Recently the idea of QND measure- 
ments has been also applied to solid-state mesoscopic 
structures (see, e.g., [7,8]). Among other recent devel- 
opments (total number of papers on QND measurements 
is over half a thousand) let us mention the experiment on 
atomic spin-squeezing using the QND measurement and 



real-time quantum feedback [9] . In this paper we discuss 
squeezing of a nanomechanical resonator using the QND 
measurement and quantum feedback. 

Recent advances in fabrication of nanomechanical res- 
onators [10-14] make possible the direct observation of 
their quantum behavior in the nearest future. For the 
resonator frequency loqI^-k exceeding 1 GHz [12] the con- 
dition T < fiLOQ (we use fcs = 1) is satisfied at temper- 
ature T below ~ 50 mK. Even in the case T ^ fiLOi^ 
the quantum behavior is in principle observable [1] when 
Tt/Q < fi/2, where Q is the resonator quality factor and 
r is the typical observation time. This condition can be 
well satisfied even for MHz-range resonators for r com- 
parable to the oscillation period, i.e. if we can monitor 
the oscillations with the measurement bandwidth on the 
order of loq, as in Ref. [14]. There is a rapid experimen- 
tal progress in monitoring the oscillating position of a 
nanoresonator using radio- frequency single-electron tran- 
sistor (RF-SET) [13,14] or quantum point contact (QPC) 
[15] (at present RF-SET seems to be much more effi- 
cient). In particular, the position measurement accuracy 
Aa; within the factor 5.8 from the standard quantum limit 
Aa;o has been demonstrated [14]; here Axq = y/fi/2rnujQ 
is the width (standard deviation) of the ground state of 
the oscillator with mass m. Measurement of the nanores- 
onator position by RF-SET or QPC has also received a 
significant theoretical attention [16-22]. The process of 
measurement transfers the energy from the detector to 
the nanoresonator leading to its "heating" [16,17]. A 
possible way to prevent such heating is using the quan- 
tum feedback control of the nanoresonator [22] (another 
idea for cooling has been proposed in Ref. [23]). 

The quantum feedback of mesoscopic solid-state sys- 
tems is a relatively new subject [24] , though in quantum 
optics the quantum feedback has been proposed more 
than a decade ago [25] and has been already realized 
experimentally [9]. The feedback analyzed in Ref. [22] 
assumes continuous monitoring of the nanoresonator po- 
sition with constant "strength" of measurement and al- 
lows cooling of the nanoresonator practically down to 
the ground state. However, it does not allow squeez- 
ing of the nanoresonator state (below Aa^o), that would 
be important, for example, for the ultrasensitive mea- 
surement of a force acting on the nanoresonator. (More 
correctly, squeezing in Ref. [22] is possible only in the un- 
realistic case of a strong coupling between the nanores- 
onator and detector.) In this paper we analyze the case 
of the nanoresonator monitoring with the measurement 



strength modulated in time (for example, modulating the 
bias voltage of the QPC or RF-SET), basically applying 
the idea of stroboscopic QND measurements [2,3] to the 
nanoresonator. We show that significant squeezing of the 
nanoresonator state can be achieved when the modula- 
tion frequency w is close to 2ujQ/n, n = 1, 2, . . ., similar 
to the usual QND result. The difference from the case 
of Refs. [2,3] though is that we consider continuous in 
time measurement and assume weak coupling between 
the nanoresonator and the detector; also, we use the 
quantum feedback to prevent the heating due to mea- 
surement. We would also like to notice the difference 
between our proposal and squeezing of a nanoresonator 
proposed in Ref. [26], which is a scaled-down version of 
the classical thermomechanical noise squeezing [27] using 
the modulation of the resonator spring constant. 

II. SYSTEM AND MODEL 

For simplicity we consider the nanoresonator measured 
by the low-transparency QPC (though our results are ap- 
plicable to the RF-SET as well) and the system Hamil- 
tonian is 



H = Ho + Ti-de 



'Hi: 



Hfb, 



(1) 



where the first term describes the oscillator (nanores- 
onator) : 



,2;P,2 



Ho 



2m 



(2) 



(p and X being the momentum and position operators), 
the last term 



H 



fb 



-Tx 



(3) 



describes the feedback control of the nanoresonator by 
applying the force T(t\ while H.det and H.int correspond 
to the detector and its interaction with the nanoresonator 
similar to Refs. [28,17]: 



+ ^ (Ma|ar + H.c.), 



(4) 
(5) 



Here a^ ^ and ai^r are the creation and annihilation op- 
erators for two electrodes of the QPC, and for simplic- 
ity we assume no relative phase between the tunneling 
amplitudes M and AM (taking this phase into account 
is simple [29-31], but makes the formalism significantly 
lengthier). For a given position state \x) of the os- 
cillator, the average detector current is Ix — 2Tr\M + 
AMxl'^piPre'^V/h, where V is the QPC voltage which 



may vary in time with frequency lo comparable to Wq, e 
is the electron charge, and pi^r are the densities of states 
in the electrodes. 

We assume a weak response of the detector, \Ix—Ix' \ ^ 
Ix + Ix' , and therefore the linear dependence of the de- 
tector current on the measured position 



/q + kx, 



(6) 



also we neglect the dependence on x of the detector cur- 
rent spectral density Sj which is assumed to be flat in the 
frequency range of interest. Because the voltage V varies 
in time, Iq, k, Ix, and Sj also depend on time, which will 
be taken into account explicitly in the next Section. 

Somewhat similar to the case of qubit measurement 
[32], we define the dimensionless (time-dependent) cou- 
pling as 



C = 



hk^ 



SimuiQ 



2 ' 



(7) 



which can also be expressed as C = A/ujoTrm where Tm — 
2Si/{kAxo)'^ is the "measurement" time which would be 
necessary to distinguish (with signal-to- noise ratio of 1) 
two position states separated by the ground state width 
Axq = ^h/2inujQ. We will mainly consider the case of 
weak coupling, C <C 1, which corresponds to a realistic 
experimental situation. As an example, C is on the order 
of 10^'^ for the parameters of experiment [14]. 

To describe the dynamics of the quantum measurement 
process, we apply the Bayesian approach [32], which is 
practically equivalent to the approach of quantum trajec- 
tories used, e.g., in Refs. [22,25,30,33]. Therefore we need 
to use the usual assumptions of the Bayesian approach; 
in particular, we assume that the internal dynamics of 
the detector is much faster than the oscillator dynamics 
(this requires eV ^ hcoo), and we assume quasicontinu- 
ous detector current (which requires lo/e 3> loq and even 
stronger inequality kAxg/e '^ luq). 

Applying the Bayesian approach to our system, we de- 
rive (derivation will be presented elsewhere) the following 
equation for the evolution of oscillator density matrix p 
in x-basis (in Stratonovich form): 

—i 1 

p{x,x') = -^[Hq +Hfb,p]x.x' + p{x,x')—- 
n bi 



x\l{t){Ix+Ix'-2{I)) 



t2 I 72 
^x I -'rr' _ / t2\ 



, (8) 



where the first term is a usual evolution due to Ho + Hfb, 
while the second term describes the evolution due to mea- 
surement, and therefore depends on the noisy detector 
current I(t); we introduced notations (/) = J Ixp{x, x)dx 
and {P) — J I^p{x,x)dx. Notice that Eq. (8) actually 
does not require the current linearity (6). Also notice 
that Eq. (8) coincides with the similar equation for the 
case of arbitrary number of entangled qubits measured 



by an ideal detector [31], if x is replaced by the index 
corresponding to the state of qubits. 

Equation (8) allows us to monitor the oscillator density 
matrix p using the measurement record I{t), while for 
simulations /(t) may be replaced with 



m = {i)+m, 



(9) 



where ^(t) is a white noise with spectral density Si. 

Translating Eq. (8) from Stratonovich into Ito form, 
using current linearity (6), and taking into account quan- 
tum efficiency 77 of the detector [32] (QPC is an ideal de- 
tector, 77 = 1, while 77 < 1 can be used for the RF-SET 
as a detector), we obtain 



p{x,x') 






+ — {x + x'-2{x))p{x,x')^{t), 



(10) 



where (x) — J x p{x,x)dx. Equation (10) is similar to 
equations derived in many publications (e.g., in Refs. 
[34,35,33,22]) for measurement of a mechanical oscilla- 
tor. Averaging Eq. (10) over the measurement record I{t) 
eliminates the last term of Eq. (10) [in Ito form averaging 
over the noise ^(i) is equivalent to using ^{t) — 0], and 
leads to the ensemble averaged equation derived in even 
larger number of papers, including, e.g., Rcf. [17]. Notice 
that the second (dccohcrcncc) term in Eq. (10) can also 
be rewritten in a standard double-commutator form (see, 
e.g., Refs. [36-38,35,39,17,33,22]) since ix~x')^p{x,x') = 
[x, [x,p]\x,x'- The nanoresonator evolution described by 
Eq. (10) does not depend on the environment tempera- 
ture because we essentially assume large (infinite) quality 
factor of the nanoresonator, so that the interaction with 
the thermal bath is much weaker than interaction with 
the detector which has infinite by assumption effective 
temperature Td - eV/2 [17,16]. 



III. QND SQUEEZING OF THE 
NANORESONATOR 

A. Modulation of the measurement strength 

We assume periodic modulation of the voltage across 
the QPC detector, V ~ f(t)Vo, which leads to the corre- 
sponding modulation of the parameters in Eqs. (8) and 
(10): 

k^f{t)ko, Ix = f{t)iIoo + kox), Si = \f{t)\So, (11) 

while quantum efficiency rj is assumed constant (in gen- 
eral case /(t) may become negative). Notice that the 
noise S,(t) has implicit time dependence because of mod- 
ulated in time spectral density Sj. The dimensionless 
couphng is modulated as C = \f{t)\Co- 



In this paper we will consider two types of modula- 
tion with frequency lo: harmonic modulation with 100% 
depth 



f{t) = (1 +cosujt)/2 



(12) 



and the square-wave (stroboscopic) modulation with 
pulse width St 

^(,)^ri, \t-j.2nH<St/2,j^l,2,... 

■^ ^ ' \0, otherwise . ^ ' 

Notice that |/(t)| < 1, so fco and Co correspond to the 
maximum coupling. Since f{t) reaches zero in both types 
of modulation, the conditions eV 3> hujQ and kAxo/e 3> 
luq required for the Bayesian formalism are violated dur- 
ing a fraction of the modulation period. However, the 
Bayesian formalism is trivially correct at V — when the 
nanoresonator is not measured (neglecting small remain- 
ing coupling [17]). So, assuming that the intermediate 
regime for which the Bayesian formalism is inapplicable, 
is realized during only a short fraction of the modulation 
period, we still use Eq. (10) for the analysis. 



B. Simplified equations for the Gaussian states 

Following Refs. [40,41,33,22], we assume that the os- 
cillator state is Gaussian [42]: 



p{x,x') 



X exp 



X exp 



V2¥A 



exp 



2Dx 



{x - x'f {D^D^ - D%) 



8D, 



h^/4: 



,(.-..)(M+(£±£:-(.,)^ 



(14) 



and therefore is characterized by only five parameters: 
average position (x) = (x) and momentum (p) = (p), 
their variances D^ = {xP') — (.t)^ and Dp = {p^) — (p)^, 
and the correlation D^p = {xp + px)/2 — {x){p). These 
parameters satisfy the generalized Heisenberg inequality 
[43] 



D^Dr, 



Dip > h'/A, 



(15) 



which reaches the lower bound for the pure states. The 
assumption of the Gaussian state can be justified by the 
fact that a Gaussian state remains Gaussian in the pro- 
cess of continuous measurement [41] (we have checked 
this statement for nonideal detectors including "asym- 
metric" detectors and for varying in time strength of mea- 
surement) and by the fact that the thermal state (natural 
initial condition) is Gaussian [42]. 

For Gaussian states Eq. (10) transforms into equations 






(16) 



(p) = -ma;2(x) + ^ sgn[/(t)] D^p ^t) + T , (17) 

Do 



K = -D^,-^-^\!{i)\Dl. 
vn Dq 



(18) 



b^ = -2^ujlD^,^^\!{t)\^^\!{t)\Dl^, (19) 
^Do?7 Do 

1 2fc2 

i>xp = -^p - mcj^i?, - — ^ |/(t)| i?,i?,p , (20) 

which practically coincide with the equations derived in 
Refs. [40,33,22], except for the time dependent j(t). It 
is interesting to notice that while Eq. (10) is a nonlinear 
stochastic equation, for which the Stratonovich and ltd 
forms are significantly different, there is no such differ- 
ence for Eqs. (16)-(20), so they can be treated as simple 
ordinary differential equations. 

Notice that the equations for Dj,, Dp, and D^^ do 
not depend on noise ^(t) and feedback force T ■, and are 
decoupled from the remaining equations. Therefore the 
evolution of the "wave packet width" a/D^ is determin- 
istic. Analyzing the possibility to squeeze the nanores- 
onator state, we will consider separately the squeezing 
of the packet width and the contribution D{x) to the to- 
tal position variance due to fluctuating position of the 
packet center (x). As will be discussed below, D^ may 
be made significantly smaller than the ground state vari- 
ance Axg using modulation / (t), while D{x) can be made 
even smaller using the feedback. 



C. Packet width squeezing 

Let us use the natural normalization of D^ and Dp by 
the ground state parameters, dx = D^/ {fi/2mujQ), dp = 
Dp/{hmuJo/2), and similarly dxp = Dxp/{h/2). Then 
Eqs. (18)-(20) can be rewritten as 



dx/cJo^2dxp-Co\f{t)\dl, 


(21) 


dp/cuo = -2dxp + (Co/ry) |/(t)| - Cq |/(t)| d% , 


(22) 


dxp/ujQ ^ dp - dx - Co \f{t)\ dxdxp ■ 


(23) 



We have analyzed these equations numerically for the 
harmonic (12) and stroboscopic (13) modulation f{t) for 
several values of the maximum coupling Co, concentrat- 
ing on the range Cq < 1. Notice that for the strobo- 
scopic modulation the evolution during each period of 
modulation can be calculated analytically using Riccati 
equations [33] that significantly simplifies the numeri- 
cal calculations. As anticipated, we have found that 
irrespectively of the initial conditions, Eqs. (21)-(23) 
approach the asymptotic solutions which oscillate with 
the modulation frequency lo. Even for small coupling, 
Co <IC 1, the asymptotic oscillations are significant in 
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FIG. 1. Dependence of the packet width squeezing 5 (max- 
imized over the modulation period) on the frequency u; of the 
harmonic modulation (12) of the measurement strength, for 
three values of the quantum efficiency rj of the detector. Solid 
lines show the numerical results, dashed lines are the analyt- 
ical results corresponding to Eqs. (30) and (32), and dotted 
lines are the asymptotes <S — y/7j. 



the case of resonance: lo ~ 2a;o for harmonic modula- 
tion and u ~ 2LUo/n for the stroboscopic modulation 
(notice that at Cq = the variances oscillate with fre- 
quency 2ciJo)- During the oscillation period the asymp- 
totic solution for dx{t) reaches the values both above 
and below the stationary solution for /(t) = 1 which is 
[33,22] dx = (\/2/Co)[(l +CoVry)i/2 - l]i/2 and becomes 
dx = 1/\A7 for Co <C 1. Most importantly, the squeezed 
state, da; < 1, may be achieved for both harmonic and 
stroboscopic modulation (momentum squeezing is also 
achieved; however, we do not analyze it in this paper). 

Figure 1 shows the maximum squeezing over the 
oscillation period for the asymptotic solution, S = 
niaxt[l / dx{t)] = max([Axo/-Dj,(i)], as a function of the 
modulation frequency u for the harmonic modulation 
(12) in the case of weak coupling Co = 0.1. One can see 
that for the ideal detector, ry = 1, the squeezing S « 1.73 
is achieved at lo = 2uj() and decreases to 5 w 1 (which 
corresponds to the ground state width) away from the 
resonance. The resonances at lo ~ 2L0()/n, n > 2, are 
barely visible and lead to small shoulders rather than 
to peaks. For nonideal detectors, ?] < 1, the height of 
the resonance peak decreases, S{2loo) ~ 1.73^^7^, while 
the width increases; the squeezing becomes impossible, 
5 < 1, at ry < 1/3. 

Much stronger squeezing of the packet width can be 
achieved for the stroboscopic modulation (13). (Corre- 
spondingly, the efficiency rj can also be much lower - see 
analytics below.) Figure 2 shows S{uj) for the ideal de- 
tector with Co = 0.5 and pulse duration 6t = O.OSToj 
where To = 2'k Jloq is the nanoresonator period. One can 
see that the sharp resonances at w = 2ujQ/n have equal 
height; however, their width decreases with n. [If mod- 
ulation (13) is modified so that f{t) = const > during 
"off" phase, then the peak height also decreases with n] 
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FIG. 2. Numerical results for the packet width squeezing 
<S as a function of modulation frequency to for the strobo- 
scopic measurement modulation (13) with the pulse width St. 
Efficient squeezing occurs at w « 2u>o/n. The height of the 
squeezing peaks is proportional to ^{5t)~^ [Eq. (41)] while 
their width is proportional to Co{5ty*/n^ ^ [Eq. (42)]. 



For smaller coupling Cq the peak height remains practi- 
cally the same, but the peak width decreases; this is the 
reason why we chose relatively large coupling in Fig. 2 
in order to have a noticeable peak width. For smaller 
pulse duration St, the squeezing peaks in Fig. 2 would 
become higher and narrower, while the detector nonide- 
ality makes peaks lower and wider (this will be evident 
from the analytical results discussed below). 



1. Evolution of the state purity 

Before discussing the analytical results for squeezing, 
let us briefly discuss the evolution of the state purity, 
Tr(p2) 



xp 



1 



where u 



O'^r.O'ri 



(P 



{n/2)/^D,D.^ 
From Eqs. (21)-(23) it is easy to derive the 



equation 



u^uJnCo\f{t)\dx{'q ^-u). 



(24) 



Since Co and d^ are positive, the asymptotic solution of 
this equation is obviously u — 1/rj and therefore the state 
purity reaches the asymptote Tr(/92) = y/r]. In particular, 
in the case of ideal detector, rj — \, the state eventually 
becomes pure (similar to the case of a qubit measurement 
[32]). It is interesting to note that the typical purification 
time is comparable to the time of reaching the asymptotic 
regime. 



2. Analytics for harmonic modulation 

Without measurement, Co — 0, Eqs. (21)-(23) have 
the solution 



v). 


(25) 


v), 


(26) 




(27) 



dx{t) = Vry-i + A2 - A cos{2ujot 
dp{t) = \/r)-^ + A^ + A cos{2ujot - 
dxp{t) ~ A sin(2ciJo^ + '/') i 

with arbitrary amplitude A and phase </?. (Notice that 
these equations satisfy the condition u = l/r).) For weak 
coupling, Co/rj <C 1, and harmonic modulation (12) in 
the vicinity of the resonance, lu ~ 2a;o, it is natural to 
look for the asymptotic solution of Eqs. (21)-(23) in the 
form (25)-(27) with 2a;o replaced with lu (actually, A 
and If vary in time with frequency lu, but variations are 
neghgible at Co/?7 ^ 1). 

To find A and tp, we substitute these equations into the 
equation JZ^'I^ f{t){r]^^ — d^ — d^ ) dt = which follows 

from the stationarity condition, /_ ; (d^ + dp) dt = 0, 
and Eqs. (21)-(22). This gives us the relation 



A^\^ 



r]- 



A^ cos If . 



(28) 



We find numerically that (p = at the resonance, lo = 
2luo. (This is quite natural, corresponding to smaller d^ 
at larger measurement strength, and is also proven be- 
low). Then from Eq. (28) we find A — l/y/3rj and there- 
fore 



Si2LUo) 



3?7 



(29) 



since the maximum squeezing S and the amplitude A are 
related as 



S = r]iA+^A^+7^-^). 



(30) 



This result confirms the numerical result for the peak 
height in Fig. 1. 

To find the shape of the resonant peak, we need one 
more equation relating A and (p. It can be obtained 
by deriving equation for dxp{t) from Eqs. (21)-(23), and 
equating the sin(u;i -I- f) component for the two sides of 
the equation (assuming Cq/t] <C 1 and w ~ 2a'o). In this 
way we obtain 



(4ciJq — uj ) a = rj Co ciJq sin f . 



(31) 



In particular, this proves that ip — Q dX lu ^ 2lU{). Com- 
bining Eqs. (28) and (31) we find the amplitude A as 



A{uu) = 



2M 



3 + g{uj) + vrH+WM + 9 



(32) 



where g{uj) — 16r7(2 — 0^/0)0)^/^0. This result gives us the 
analytical expression for squeezing S via Eq. (30). The 
corresponding squeezing is shown by the dashes lines in 
Fig. 1 , which practically coincide with the solid lines rep- 
resenting the numerical results. Notice that the lincwidth 
of the peak is proportional to CqJ ^fr\; away from the reso- 
nance A decreases to zero, and the squeezing approaches 



S = Y^, which is the same as for the case without mod- 
ulation [33]. The analytical result for S{u!) works well 
for coupling Co up to approximately 0.3. It is curious 
that rather complex shape of the resonance peak given 
by Eqs. (30) and (32) is quite close to the square root of 
the Lorentzian shape: 



S{lu) « V^ 1 + 



V3-1 



^ ^l + [{uj-2ujo)/Auj]y 
with Auj ~ 0.36 o^oCq/y^. 

3. Analytics for stroboscopic modulation 



(33) 



In the case of stroboscopic modulation (13) of the mea- 
surement strength, the variances dx, dp and dxp should 
follow Eqs. (25)-(27) during the "off" phase of the mod- 
ulation, while during the measurement pulse of duration 
dt ("on" phase) the parameters A and (p slowly change 
(we again assume the weak coupling limit) in accordance 
with Eqs. (21)-(23). In particular, close to the nth res- 
onant peak of Fig. 2, w « 2uJo/n, the phase ip should 
change during the pulse by the small amount 

Sip = — 2ti>o \- 2Trn Ri Tin {uj/ujo — 2/n) (34) 

UJ 

in order to match 2Tr/uj periodicity of the asymptotic 
solution with the periodicity of free oscillations (25)-(27). 
On the other hand, S(p can be found from the equation 



ip = -4ujoCoV''\f{t)\dxp/[{dp - dxf + Ml 



(35) 



which follows from from Eqs. (21)-(23). Integrating Eq. 
(35) within the pulse interval |i| < dt/2 using Eqs. (25)- 
(27) in which A and ip are assumed constant, we obtain 
Sip = —Cosin(uj()St)/riA. Combining this result with Eq. 
(34) we obtain an equation relating A and ip: 

nn A{uj/ujo — 2/n) — 77^ Co sin(a;o'5t) simp. (36) 

To obtain one more equation for A and p, we use the 
condition J_sf/2{dx + dp) dt — 0. Expressing the deriva- 
tive dx + dp from Eqs. (21)-(22) and using Eqs. (25)-(27) 
we get the equation 



AbJaSt ^ y^rj ^ + A'^ sin(u'o'^i) cos 93. 



(37) 



Equations (36) and (37) are sufficient to find A for the 
nth resonance, though the expression is quite long: 



AHco) 



2rj ^ sin (uioSt) 



(uj) + Jb^{uj) + Agiuj) siri^iLUoSt) 



(38) 



B 



where B{uj) = g{uj) + {dJoSt)"^ — sin {ojnSt) and (/{oj) = 
7r^n^(2/n — a;/a;o)^r?/Co- The squeezing S is obtained 
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FIG. 3. Squeezing for the stroboscopic modulation for 
three values of the coupling with detector Co. Solid lines 
show numerical results, dashed lines (practically coinciding 
with the solid lines) are the analytical results given by Eqs. 
(38) and (30), and the dotted lines are calculated using the 
simplified equation (40). 



from this result using Eq. (30). The corresponding ana- 
lytical curves are plotted in Fig. 3 by the dashed lines 
which practically coincide with the numerical results 
shown by the solid lines. One can sec that the analytics 
works well even for Co = 1, even though it was assumed 
Co <C 1 for the derivation. 

The value of squeezing at lu — 2LUo/n (peak height) can 
be obtained from Eq. (38), but it is easier to use Eq. (37) 
with ip = [which follows from Eq. (36)], that leads to 
the result 



S{2uJo/n) = y/rj. 



I ujoSt + sm{LLioSt) 
ujoSt — sm{ujoSt) 



(39) 



The analytical results significantly simplify in the case 
of short pulses, St ^ Tq = 27r/a;o, then 



A'iu;) = 



6/{ioo6t)^T] 



1 



1 



67r y/rj n^ (lo — 2ujQ / n) 



(40) 



which corresponds to the peak squeezing 

5(2(jJo/n) = 2y^/LL!o5t 
and the full width at half height of S{uj) equal to 

Aw = 'iCaiSt)^cJn/'^n^V^- 



(41) 



(42) 



The curves calculated using equation (40) arc shown in 
Fig. 3 by the dotted lines. As one can see, there is a no- 
ticeable difference from the numerical results away from 
the resonance; however, the main part of the peak is fit- 
ted quite well. 

Equation (41) shows that the maximum squeezing in 
our model does not depend on coupling Co. Nanorcs- 
onator interaction with extra environment (for example, 
due to finite quality factor) would obviously change this 



conclusion, because the corresponding decay of squeez- 
ing would compete with the squeezing "build up" , which 
rate is propotional to Co as follows from Eqs. (21)-(23). 



D. Quantum feedback of the packet center 

While the width of the monitored Gaussian packet can 
be squeezed below the ground state width as shown in 
the previous subsection, the center of the packet under- 
goes random evolution described by Eqs. (16)~(17), and 
without feedback diffuses far away from the origin. For 
our model which assumes infinite quality factor of the 
nanoresonator, the packet center evolves infinitely far 
away because the back-action from the detector "heats 
up" the nanoresonator to formally infinite effective tem- 
perature (voltage) of the detector [33,17,22]. To prevent 
deviation of the packet center from the origin, we can 
apply the quantum feedback described by the force JF in 
Eq. (17). Similar to Refs. [33,22], we choose it as 



T = -mujojx{x) - jp{p). 



(43) 



Notice that the random oscillating evolution of (x) and 
(p) can be extracted from the measurement record using 
Eqs. (16)~(17) even when the measurement is performed 
during only small fraction of the period. 

Following Refs. [33,22], we characterize the distribu- 
tion of the packet center position (x) and center mo- 
mentum (p) by the ensemble averages (over realizations) 
((x)) and {{p)) and the variances D(^x) — ii^)'^) ~ {{^))'^i 
Di,) - {{P?) ((P»^ i?(.)(rt - ((x)(p)) - ((x))((p)). In 
the notation of doubled angle brackets the inner brackets 
mean averaging with individual density matrix p, while 
the outer brackets is averaging over realizations. 

The equations for ((x)) and ((p)) derived from Eqs. 
(16)-(17) lead to the ensemble averaged position evolu- 
tion as 



((i» + 7p((i» + (^0 + lxOJa){{x)) - 0, 



(44) 



from which it is clear that ((x)) relaxes to zero for positive 

Ip- ^ 

Introducing dimensionlcss 

variances di^^) = D(^x^2mujQ/h, d(^p^ = D(^p'^2/hmu}o, and 
di^x){p) = D(^x^p'^2/h, we derive the following equations 
from Eqs. (16)-(17): 

d{x)/t^o = 2d(j.)(p) +CQ\f{t)\dl, (45) 

d{p)/u}Q = -2d(j.)(p) - 2^F(i(j.)(p) - 2Fd(p) 

+Co|/(t)K, (46) 

d{x){p)/^o — d{p) — d/^x-j — iiFdi^x) — Fd(^x)(p) 

+Co|/(t)|d,d,p, (47) 

where F = 7p/wo and /i — ^x/lp are the dimensionlcss 
feedback parameters. 



We have simulated these equations numerically using 
the asymptotic solutions of Eqs. (21)-(23) for dx, dp, and 
dxp- We have mostly studied the resonance lo — 2ujq in 
the weakly coupling regime. The main finding is that for 
both harmonic (12) and stroboscopic (13) modulation 
of measurement, the center position variance du) can 
be made much smaller than the packet variance dx at 
time moments t = j27r/aj when the packet squeezing is 
at maximum. Therefore, the corresponding worsening of 
the ensemble-averaged squeezing defined as [dx + di^x) + 
((x))^/Axg)~^ is negligible. 

As an example, for the stroboscopic modulation with 
6t/To — 0.05 and Co = 0.1 at cj = 2a;o, the ratio d(x)ldx 
at t = jT'o/2 is around 0.4% for ,u = and F = 0.5, and 
around 0.1% for /i = 10 and F 3> 1. [The term with ^x in 
Eq. (43) is not necessary; however, it improves squeezing 
of the packet center, so nonzero /x is beneficial.] The ratio 
d{x) /dx decreases with decrease of the pulse width 5t and 
decrease of coupling Cq. It is important to notice that 
d(^x) scales linearly with Co as follows from Eqs. (45)- 
(47), while dx as well as dp and dxp do not depend on 
Co at w = 2a;o and Co/?? <C 1 (this was also checked 
numerically). Therefore the ratio dix)ldx can be made 
arbitrary small at small coupling. 

The analytical results (which will be described in more 
detail elsewhere) show that in the case [iioobt 3> 1 and 
F 3> 1 , the variance of the packet center at the middle of 
the measurement pulse is dix\ = Co(tiJo(5t)^/24^?7, which 
can be obviously made much smaller than S~^ given by 
Eq. (41) at sufficiently small Co/^/rj, uJodt, and/or suffi- 
ciently large /i. 



E. Observability of the squeezed state 

The fact that the squeezed state of a nanoresonator can 
be prepared by the modulated measurement and quan- 
tum feedback, does not automatically mean that this 
state may be useful for the measurement of extremely 
weak forces, and even that such state can be checked ex- 
perimentally in a straightforward way. As an example of 
such problem, in one of setups analyzed in Rcf. [44] the 
squeezed in-loop optical state is realized by using quan- 
tum feedback, but the squeezing of the output light is 
impossible. 

We have studied the possibility to verify the squeezed 
state of the nanoresonator in the following way. After the 
preparation of the squeezed state by stroboscopic mea- 
surement and feedback, the feedback at some moment 
{t = 0) is switched off, while the stroboscopic measure- 
ment continues. Considering for simplicity the case of one 
measurement per nanoresonator period (n = 2, w = ujq), 
we calculate the average of the position measurements 
(each pulse gives a very imprecise measurement because 
of weak coupling) : 



N 



]To+St/2 



N f-{Stko JjTo-5t/2 



(48) 



The idea is that for a squeezed initial state, Xn can 
be much closer to zero than if we would start with the 
ground state. 

The analysis of the distribution of Xj^ (over realiza- 
tions) is very simple in the case of instantaneous but 
imprecise measurements, St -^ 0, Co5t — const, since the 
Hamiltonian evolution of the resonator in between the 
measurements can be completely neglected, and therefore 
N measurements are equivalent to one iV-times stronger 
measurement. This gives us the variance of Xn equal to 



D 



X,N 



1 



1 



2muJo \S NCqujoSI 



(49) 



where the first term is due to the initial packet width, 
while the second term is the inaccuracy of the mea- 
surement which improves with N . Obviously, at N ^ 
l/CiiLOiiSt this variance for a squeezed state {S > 1) is 
significantly different from the variance for the ground 
state {S ~ 1). Even though this difference can be rigor- 
ously verified only by performing many experiments to 
accumulate statistics for Dx.Ni it can be observed even 
in the single experiment with good reliability if 5 3> 1. 
(In a single realization the failure probability for distin- 
guishing squeezed and ground states is crudely S^^'"^.) 

Unfortunately, this result requires the assumption of 
infinitely strong coupling, so it is not obvious if it holds 
in the case of weak coupling or not. The possible prob- 
lem is that for sufficiently large N which makes the 
second term in Eq. (49) sufficiently small, the nanores- 
onator heating due to measurement back-action may al- 
ready eliminate the squeezing (the feedback is off). We 
have calculated Dx.n for stroboscopic modulation nu- 
merically using Eqs. (16)-(17) and have found that there 
is still a range of N where the squeezed and ground ini- 
tial states lead to significantly different Dx,n and there- 
fore can be reliably distinguished. As an example, for 
Co = 0.1, rj — 1, and 5t — 0.02To, the normalized variance 
dx,N = Dx,N'2mLUo/h achieves a minimum of 0.078 (at 
TV ~ 4 X 10"^), which is crudely two times larger than the 
contribution from the initial squeezing 1/5 = 0.036, and 
is still significantly smaller than the ground state limit 
dx,N > 1- We have found numerically that the mini- 
mum dx.N scales linearly with the pulse width 6t (simi- 
larly to 1/5) and practically does not depend on coupling 
Co at Co < 1. This hints that the product S x imnxdx.N 
is practically a constant approximately equal to 2, and 
therefore verification of the squeezed state by a weakly 
coupled detector is almost as efficient as the verification 
by instantaneous measurements (within a factor of about 
2). 



IV. CONCLUSION 

In this paper we have shown that the uncertainty of the 
nanoresonator position can be squeezed significantly be- 
low the ground state level by the modulated in time mea- 
surement of the nanoresonator position with the QPC or 
RF-SET detector. The measurement strength is mod- 
ulated by applying the periodic voltage across the de- 
tector. The mechanism of squeezing is similar to the 
QND measurements [1] , though a significant difference in 
our case is the continuous measurement with weak cou- 
pling to the detector. We have considered harmonic (12) 
and stroboscopic (13) modulations and found that only a 
moderate squeezing S < \/5ri is possible for the harmonic 
modulation with frequency w w 2^0 [see Eqs. (29) and 
(32), and Fig. 1]. However, the stroboscopic modulation 
can lead to an arbitrary strong squeezing S < 2y/3ri/LUo5t 
for sufficiently short measurement pulses dt applied with 
frequency w — 2ujo/n [see Eqs. (39)-(42) and Fig. 2]. 
Obviously, the state width oscillates with time, so that 
the maximum position squeezing is achieved periodically 
(with period close to tt/loo), while the maximum squeez- 
ing of momentum happens with 7:/2luo = T^/A time shift 
(when 5 < 1). 

While the modulated measurement squeezes the width 
of the state (packet), the position of the packet center {x) 
fluctuates due to random back-action from the detector; 
so to keep the packet center near a; = we need to apply 
quantum feedback. We have found that the feedback 
can keep the deviation of {x) from zero much smaller 
than the packet width, which means that the ensemble- 
averaged squeezing practically does not differ from the 
packet width squeezing. 

In this paper we have used the Bayesian formalism [32] 
for the description of the quantum measurement process. 
However, since the obtained equations practically coin- 
cide with the equations used in Refs. [33,22], we have 
followed those papers to a large extent, especially for the 
analysis of the evolution of the Gaussian states. 

An important issue is the possibility to use the 
squeezed states of the nanoresonator for the measurement 
of weak forces with the accuracy beyond the standard 
quantum limit. Even though we did not consider this 
question explicitly, we have found that the state squeez- 
ing can be verifled (with high reliability) even in a single 
measurement run by a weakly coupled detector. 

The main drawback of the present theory is the as- 
sumption of a very large quality factor Q of the nanores- 
onator. Crude preliminary analysis indicates that our 
results for the stroboscopic modulation are valid for 
Q 3> S'^/Cq^ ~ ^/Co{ujQ5tY and sufficiently small 
temperature of the environment, T/huJo <C y/rjQCa/S^ ^ 
QCo{uJoSt)^/ri. These conditions seem to be within 
present-day experimental reach for moderate squeezing. 



ACKNOWLEDGMENT 

The authors would hkc to thank D. Avcrin, A. Do- 
herty, S. Habib, K. Jacobs, K. Likharev, I. Martin, and 
G. Milburn for fruitful discussions and remarks. This 
work was supported by NSA and ARDA under the ARO 
grant DAAD 19-01-1-0491 (R.R. and A.N.K.) and by the 
National Security Agency (K.S.). 



[1] V. B. Braginsky and F. Ya. Khalili, Quantum measure- 
ment, (Cambridge Univ. Press, Cambridge, 1992). 
[2] V. B. Braginsky, Yu. I. Vorontsov, and F. Ya. Khalili, 

JETP Lett. 27, 276 (1978). 
[3] K. S. Thorne, R. W. P. Drever, C. M. Caves, M. Zim- 
mermann, and V. D. Sandberg, Phys. Rev. Lett. 40, 667 
(1978). 
[4] C. M. Caves, K. S. Thorne, R. W. P. Drever, V. D. Sand- 
berg, and M. Zimmermann, Rev. Mod. Phys. 52, 341 
(1980). 
[5] H. J. Kimble, Y. Levin, A. B. Matsko, K. S. Thorne, and 

S. P. Vyatchanin, Phys. Rev. D 65, 022002 (2002). 
[6] V. B. Braginskii, Physics-Uspekhi 46, 81 (2003). 
[7] D. V. Averin, Phys. Rev. Lett. 88, 207901 (2002). 
[8] L. Bulaevskii, M. Hruska, A. Shnirman, D. Smith, and 

Y. Makhlin, Phys. Rev. Lett. 92, 177001 (2004). 
[9] JM Geremia, J. K. Stockton, and H. Mabuchi, Science 
304, 270 (2004). 

[10] A. N. Cleland and M. L. Roukes, Appl. Phys. Lett. 69, 
2653 (1996); A. N. Cleland and M. L. Roukes, Nature 
392, 160 (1998); E. Buks and M. L. Roukes, Europhys. 
Lett. 54, 220 (2001). 

[11] D. W. Carr and H. G. Craighead, J. Vacuum Sci. Tech. 
B 15, 2760 (1997); H. G. Craighead, Science 290, 1532 
(2000); M. Zalalutdinov, B. Ilic, D. Czaplewski, A. Zehn- 
der, H. G. Craighead, and J. M. Parpia, Appl. Phys. Lett. 
77, 3287, (2000). 

[12] X. Ming, H. Huang, C. A. Zorman, M. Mehregany, and 
M. L. Roukes, Nature 421, 496 (2003). 

[13] R. G. Knobel and A. N. Cleland, Nature 424, 291 (2003). 

[14] M. D. LaHaye, O. Buu, B. Camarota, and K. C. Schwab, 
Science 304, 74 (2004). 

[15] A. N. Cleland, J. S. Aldridge, D. C. DriscoU, and A. C. 
Gossard, Appl. Phys. Lett. 81, 1699 (2002). 

[16] M. Blencowe, Phys. Reports 395, 159 (2004); M. P. 
Blencowe and M. N. Wybourne, Appl. Phys. Lett. 77, 
3845 (2000); A. D. Armour, M. P. Blencowe, and Y. 
Zhang, Phys. Rev. B 69, 125313 (2004). 

[17] D. Mozyrsky and I. Martin, Phys. Rev. Lett. 89, 018301 
(2002). 

[18] A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. 
Rev. Lett. 88, 148301 (2002). 

[19] A. Y. Smirnov, L. G. Mourokh, and N. J. M. Horing, 
Phys. Rev. B 67, 115312 (2003). 

[20] K. Schwab, Appl. Phys. Lett. 80, 1276 (2002). 



[21 
[22; 
[23; 
[24; 
[25; 

[26; 

[27 
[28' 
[29; 

[30 

[31 

[32; 

[33 

[34' 
[35; 

[36 

[37 

[38; 

[39 

[40 
[41 

[42; 

[43; 

[44; 



R. Knobel and A. N. Cleland, Appl. Phys. Lett. 81, 2258 

(2002). 

A. Hopkins, K. Jacobs, S. Habib, and K. Schwab, Phys. 

Rev. B 68, 235328 (2003). 

I. Martin, A. Shnirman, L. Tian, P. ZoUer, Phys. Rev. B 

69, 125339 (2004). 

R. Ruskov and A. N. Korotkov, Phys. Rev. B 66, 041401 

(2002); A. N. Korotkov, cond- mat/0404696 (2004). 

H. M. Wiseman and G. J. Milburn, Phys. Rev. Lett. 70, 

548 (1993). 

M. P. Blencowe and M. N. Wybourne, Physica B 280, 

555 (2000). 

D. Rugar and P. Griitter, Phys. Rev. Lett. 67, 699 (1991). 

S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 

A. N. Korotkov and D. V. Averin, Phys. Rev. B 61, 

165310 (2001). 

H. S. Goan and G. J. Milburn, Phys. Rev. B 64, 235307 

(2001). 

A. N. Korotkov, Phys. Rev. B 67, 235408 (2003). 

A. N. Korotkov, Phys. Rev. B 60, 5737 (1999); A. N. 

Korotkov, Phys. Rev. B 63, 115403 (2001). 

A. C. Doherty and K. Jacobs, Phys. Rev. B 60, 2700 

(1999). 

M. B. Mensky, Physics-Uspekhi 168, 1017, (1998). 

M. J. Gagen, H. M. Wiseman, and G. J. Milburn, Phys. 

Rev. A 48, 132, (1993). 

G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

A. O. Caldeira and A. J. Leggett, Ann. Phys. (N.Y.), 

149, 374 (1983). 

C. Caves and G. J. Milburn, Phys. Rev. A 36, 5543 

(1987). 

W. H. Zurek, S. Habib, and J. P. Paz, Phys. Rev. Lett. 

70, 1187 (1993). 

J. Halliwell and A. Zoupas, Phys. Rev. B 52, 7294 (1995). 

J. K. Breslin and G. J. Milburn, Phys. Rev. A 55, 1430 

(1997). 

C. W. Gardiner, Quantum noise, (Springer, Berlin, 

1991). 

V. V. Dodonov, E. V. Kurmyshev, and V. I. Man'ko, 

Phys. Lett. A 79, 150 (1980). 

H. M. Wiseman and G. J. Milburn, Phys. Rev. A 49, 

1350 (1994). 



